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The surface quantum Hall state, magneto-electric phenomena and their connection to axion elec- 
trodynamics have been studied intensively for topological insulators. One of the obstacles for ob- 
serving such effects comes from nonzero conductivity of the bulk. To overcome this obstacle we 
propose to use an external magnetic field to suppress the conductivity of the bulk carriers. The 
magnetic field dependence of galvanomagnetic and electromagnetic responses of the whole system 
shows anomalies due to broken time-reversal symmetry of the surface quantum Hall state, which 
can be used for its detection. In particular, we find negative linear dc magnetoresistivity and a 
quadratic field dependence of the Hall angle, shifted rf cyclotron resonance, nonanalytic microwave 
transmission coefficient and saturation of the Faraday rotation angle with increasing magnetic field 
or wave frequency. 
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I. INTRODUCTION 

Unlike ordinary band insulators, semiconductors or 
semimetals, a recently identified class of materials - topo- 
logical insulators (TIs) 1 ^- - exhibit unusual conducting 
states on sample boundaries. On the surface of a three- 
dimensional TI such a state is characterized by a nodal 
spectrum with a single Dirac cone (or, in general, with 
odd number of Dirac cones). If time-reversal symme- 
try (TRS) is broken, an energy gap A is induced at the 
Dirac points, and the surface state exhibits the anoma- 
lous quantum Hall (QH) effector— Although the genera- 
tion of a sizable Dirac gap requires an effort, the surface 
quantum Hall state in TIs is of great interest because it 
gives rise to rich magneto-electric phenomenal— specific 
to axion electrodynamics^ 

There is however a serious obstacle for identifying 
the surface-related magneto-electric phenomena in three- 
dimensional TIs. It stems from dissipative bulk conduc- 
tivity which generally cannot be ignored because of the 
complex band structure of three-dimensional TI where 
the Fermi level does not necessarily lie in the bulk band 
gap or crosses both the surface and bulk states*^ - — For 
a TI film with thickness d, bulk zero-field dc conductivity 
£ and surface QH conductivity a H ~ e 2 /h, the contribu- 
tion of the surface with respect to the bulk is character- 
ized by parameter e 2 /hdY,J£ It has been shown that the 
well-resolved surface magneto-electric effects, such as the 
Kerr or Faraday rotation, require sufficiently large val- 
ues of e 2 /hdT, ^ 10 ! 11 If, however, the bulk conductivity is 
much larger than the surface one, i.e. e 2 /hdY, <C 1, is it 
still possible to resolve surface magneto-electric effects in 
TIs? In this paper we demonstrate such a possibility on 
several different examples of electrodynamic phenomena. 

We show that the surface contribution to the electrody- 
namics of TIs becomes more pronounced when the bulk 
conductivity is suppressed by an external magnetic field 
B and by finite frequency uj of an applied ac electromag- 
netic field. This can be seen from Boltzmann transport 



theory expressions for the longitudinal H L and transverse 
(Hall) S H bulk conductivities (see e.g. Refs. Il8lfl9t ): 
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[1-icjT] 2 + n 2 T 2 ' H [i - iujr} 2 + n 2 r 2 



,(1) 



where u and the cyclotron frequency VL = eB/mc are 
both assumed much smaller than the frequency A/h as- 
sociated with the surface Dirac gap A: 



, n < A/h, 



(2) 



m is the effective cyclotron mass and r is the elastic scat- 
tering time. Clearly, with increasing |fi|r and lot, the 
real parts of conductivities E L H can be made compara- 
ble with e 2 /hd, even though for the zero-field dc case 
e 2 /hdTj < 1. Under these conditions the TRS breaking 
on the TI surface leads to anomalous galvanomagnetic 
and electromagnetic responses of the whole system. In 
particular, we find (i) negative linear dc magnetoresis- 
tivity and Hall angle quadratic with B, (ii) rf cyclotron 
resonance at shifted frequency 



dE/r|<r H | + |fi| 



(3) 



(iii) nonanalytic independence of the microwave trans- 
mission coefficient and (iv) saturation of the Faraday ro- 
tation angle with increasing magnetic field or wave fre- 
quency. Below we explain in detail how these anomalies 
are related to the surface QH state and how they can be 
used for its experimental identification. 

The paper is organized as follows. In Sec. [TT] we for- 
mulate the main equations of electrodynamics of a TI 
film and discuss approximations used throughout the pa- 
per. Then we present the solutions of this electrodynamic 
problem for different physical situations: galvanomag- 
netic phenomena (Sec. IIII[) . cyclotron resonance (Sec. 
TV|) , and electromagnetic transmission and Faraday rota- 
tion effects (Sec. EJ. Finally, in Sec. EH we summarize. 
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a time-dependent electromagnetic (EM) field (oc e~ tut ) 
adiabatically. This is justified for low frequencies lj <C 
A/h which can, at the same time, be much smaller than 
the plasma frequency (see below). In particular, the sur- 
face states remain dissipationless,?' 9 i.e. the surface cur- 
rent density induced by the electric field E(r) can be 
written as 

j s (r) = [4<5(z - d/2) + a h J{z + d/2)]z x E(r), (4) 

where z is the unit vector perpendicular to the film. The 
use of the delta functions in Eq. (@| is justified if the pene- 
tration lengths of the surface states into the bulk is much 
smaller than the film thickness d. Also, in this case there 
is no magnetically induced energy gap in the interior of 
the film. Therefore, the bulk conductivity tensor con- 
tains both the dissipative (longitudinal) T, L and Hall T, H 
components which are both w-dependent and given for 
a single (lightest) carrier group by the Boltzmann trans- 
port theory expressions (JlJ. The resulting bulk current 
density is 

j B (r) = E i E(r) + E ff zxE(r). (5) 

To find the EM field inside the TI, we use the thin-film 
approximation d <C c/w, which for the upper frequency 
limit huj = A = 10 meV implies d <C ch/A « 20 /im. 
In addition, the film thickness, d, should be smaller than 
the skin penetration depth, 5: 



FIG. 1: (a) Schematic geometry of a galvanomagnetic exper- 
iment with a TI film subject to perpendicular dc magnetic 
field B and electric current I. Magnetic fields Bt,& at the 
outer top and bottom surfaces are generated by the current 
(see, text), dc longitudinal resistivity p xx (b) and Hall angle 
pyx/ Pxx (c) versus external magnetic field B and normalized 
bulk conductivity dY</a H ; p xx and B are in units of 1/S and 
mc/er, respectively. Gray regions are dielectric media. 



II. FORMULATION OF THE PROBLEM 

TRS breaking on the surface of a TI can be achieved by 
coating it with thin layers of a ferromagnetic (FM) ma- 
terial f£ i 20 ' 21 magnetized perpendicularly to the TI film 
plane by an external dc magnetic field B (see, Fig. [TJi). 
The FM magnetization acts on the electron spin, gen- 
erating an energy gap A at the Dirac point, which can 
be described by a Hamiltonian H = ( — I) s uerp + Aa z ,^ 
where a — (a x , <T y ,<j z ) are the spin Pauli matrices, and v 
and p are the velocity and momentum for top s = +1 (t) 
and bottom s = —I (b) surfaces. It is assumed that 
the Fermi level lies within the gap A such that both 
surfaces of the TI have vanishing dissipative longitudi- 
nal conductivities and nonzero quantized Hall conduc- 
tivities^ <r^ b = e 2 vt,b/h, with half-integer filling factors 
i/ t b J^2. If the variation of B is restricted by Eq. the 
surface states remain on the QH plateaus, and er^ 6 do 
not depend on |_B|<22 

We also assume that the surface states respond to 



„ / ic 2 I — iuiT c / i 

d«<S = Re\/- — = — ReWl + — , (6) 

V 4ttuj S uj p V wt 

where ui p = Airne 2 v F /p F is the plasma frequency of the 
bulk carriers with density n, Fermi velocity v F and mo- 
mentum p F . For n = 10 14 cm -3 and v F = 0.5 • 10 6 ms" 1 
the plasma frequency is lo p rs fO 13 s , yielding the lower 
bound for the skin depth ~ c/u p = 0.3 /im [see Eq. ©]. 
Note also that the plasma frequency u> p w 10 13 s _1 is of 
the order of the frequency A/H rj 1.5 • I0 13 s _1 related 
to the surface gap A = 10 meV. Therefore, in addition 
to requirement w C A/fi we have uj -C uj p . 

Under condition ([6]), the electric field E inside the 
film can be approximated by the average value (E) = 

J d ^, 2 'Edz/d. The equation for (E) is obtained by aver- 
aging the Maxwell equation V x B(r) + (ieoj / c)E(r) = 
(47r/c)[j B (r) + j s (r)] over the film thickness (e is the di- 
electric constant): 

d(E L ~ ieu/4ir) (E) + (dS H + a H )i x (E) = 

= (c/4tt) z x (B t - B 6 ), a H = (e 2 /h) (j/ t + i*). (7) 

The external EM perturbation enters via magnetic fields 
B t j, at the outer top and bottom surfaces of the TI. Wc 
specify Bt.h in each concrete situation considered below. 

III. GALVANOMAGNETIC PHENOMENA 

We begin by considering dc galvanomagnetic phenom- 
ena in the standard four-contact geometry (see, e.g. 
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IV. CYCLOTRON RESONANCE 



Rep, 




FIG. 2: Cyclotron resonance: Real part of longitudinal ac 
resistivity p xx [see, Eq. (|10[) ] versus frequency uj and mag- 
netic field B (in units of 1/r and mc/er, respectively); 
dS/a H =10. 



Ref. [TH and Fig. QJi) in the presence of perpendicular 
magnetic field B and electric current I. The current 
induces the jump of the magnetic field across the film, 
z x (B t — Bf,) = (in/cw)!, so that Eq. ([7|) determines 
the longitudinal and Hall electric fields in terms of given 
/: (E x ) — p xx I/dw and (E y ) = p yx I/dw, where the 
longitudinal p xx = E^E 2 + [E H + o H jdf) and Hall 
Pyx — — Pxx[^ H + o , H /d]/'S L resistivities are given by 



1 

Pxx g 
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dE 
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(8) 
(9) 



Let us consider now the bulk cyclotron resonance. It 
can be realized in a contactless setup where the sample 
is placed in the maximum of the electric field of an rf 
resonator normal mode, which generates an antisymmet- 
ric magnetic field across the sample, i.e. B t = B^. At 
the rf frequencies the displacement current contribution 
euj/4ir is usually smaller than conductivity E L , so that B t 
induces an average ac current density (j) — czx B t /2-7rd 
in a contactless way. The relevant observable is the lon- 
gitudinal ac resistivity which we find from Eq. (J7J as 



Pxx(u) 



(1 - iuT)/E 



1 



(10) 



For ujt > 1 it has a resonance in both u) and B dependen- 
cies, shown in Fig. [5J Note that the linear on- resonance 
relation between uj and B is the hallmark of the cyclotron 
resonance. However, the resonant frequency ([3]) is shifted 
with respect to because of the additional drift in the 
direction of the current, induced by the surface contribu- 
tion to the Hall electric field oc er H /dE in Eq. ([9]). For 
bulk Drude conductivity E = ne 2 v F T /p F the frequency 
shift depends on the bulk carrier density n, Fermi veloc- 
ity v F and momentum p F — hk F as 



.14 



2wndv F /\v t + v b \k f 
~ 3 , d = 50 nm, v„ = 



101 



\v t + 
10 12 



Vb\ 



= 1 the resonance frequency shift is u> res - 
i.e. well below both A/h w 1.5 • 10 13 s 



(11) 
_1 and 
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For a H — we recover the usual results of the magne- 
totransport theory: p xx = 1/E and (E y ) / (E x ) = —fir. 
The details of the method used in our calculations are 
very well described in literature (see, e.g. Refs. [l8llT9l ). 
The -B-field independent resistivity p xx = 1/E reflects 
strong cyclotron drift in the direction of the current (x), 
induced by the crossed magnetic B and electric Hall E y 
fields. Moreover, this conclusion remains valid in the 
nonlinear electrodynamics where the dependence of the 
bulk conductivity on the magnetic field of the current (or 
of an external EM wave) is taken into account^ 

However, according to Eq. @ the nonzero surface con- 
ductivity generates an additional Hall field oc a H /c?E that 
affects the cyclotron drift in the direction of the current. 
For this reason p xx (jHJ acquires the magnetic field de- 
pendence. Moreover, since on the QH plateau a H does 
not change with \B\, resistivity p xx (B) exhibits linear 
behavior, whereas the Hall angle has the anomalous 
quadratic B term [Figs. [TJd and c]. We note that the 
anomalous terms in p xx and p yx remain identifiable even 
if several bulk carrier groups are taken into account, be- 
cause in that case the bulk resistivities are still regular 
analytic functions of B which can be subtracted from the 
total p xx and p yx . 



V. ELECTROMAGNETIC TRANSMISSION 
AND FARADAY ROTATION 

We now turn to microwave spectroscopy which also 
allows one to probe the surface states in Tls 16 . We con- 
sider an EM wave incident normally at the bottom sur- 
face of a TI and will analyze both the transmission co- 
efficient and the Faraday rotation of the EM field plane 
in the transmitted wave [see, also Fig. This situa- 
tion involves a new conductivity scale, viz. the inverse 
impedance of the dielectric media surrounding the TI, 
Z^ 1 = c(y/et + ^/eb)/^n, where et,b are the dielectric con- 
stants of the top and bottom materials (shown in gray in 
Fig. [3^). The new conductivity scale Zq 1 is important 
because it is much larger than the surface QH conductiv- 
ity <j h : the product Zo\a H \ is a small parameter propor- 
tional to the fine structure constant a = e 2 /cfjZ&l^: 

ZoK|=2a|i/ t +i*|/(Vet + Vet)<l- ( i2 ) 

Therefore, the response of the surface state to an EM 
wave is generically rather weak. To proceed we note that 
the electric and magnetic fields on the outer surfaces of 
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the TI film are 



E b 



(E) 
E,- 



E, 



(E) 



B i = V^zxE t , (13) 
B 6 = v^zx (Ei - E r ),(14) 



where t refers to the transmitted wave on the top surface, 
whereas i and r label the incident and reflected waves on 
the bottom (b) surface. In the thin film approximation 
the electric field on each surface equals to the average 
field. Eliminating the reflected field through E r w (E) — 
Ej, we express the magnetic field difference as B t — B(, 
z x [(v^t + V^)( E > - 2^/e^Ei], insert this into Eq. 0, 
and solve it for (E). Since the incident electric field Ej 
can be regarded real, we present the solution for the real 
part of the transmitted wave: 

ReE t « Re (E) =T(E* cos6» + Ej x z sin0), (15) 

T= , J a* +aj, = arctan^-, (16) 

i + v / ^ V a ^ 

where T is the transmission coefficient, 9 is the rotation 
angle of the EM field plane with respect to the incident 
wave (Faraday angle), and a L (u>,B) and a H (u>,B) are 
real functions given by 



a, = Re 



1 + Z a d(T, L - ieu/Air) 



[1 + Z rf(S L - ieuo/Aii)] 2 + Z 2 [dT, H + a H ] 



Re 



Z (dE H +cr H/ 



[1 + Z dCE L - teuj/Ait)} 2 + Z 2 [dZ H + a H }'- 



,(17) 



.(18) 



Similar to resistivity (J5J, the TRS breaking leads to a 
nonanalytic linear B dependence of the transmission co- 
efficient (fTU)) . To illustrate this we extract the large zero- 
field value T(0) from T(B) and plot in Figs. [3]o and c 
the difference T(B) — T(0) for a H ^ (solid curves) and 
for <j h — (dashed curves). The magnetic field range in 
which T(B) — T(0) cx \B\ can be tuned by varying pa- 
rameters lot and Z^dYi. This should help in finding the 
optimal regime for observation of the predicted anoma- 
lous magnetic-field dependence of T. 

Figure [3ji shows the low-frequency Faraday angle 9 
(JTHJ) in units of the fine structure constant a as a function 
of the magnetic field B for zero and finite bulk conductiv- 
ity S. For £ = the Faraday angle contains only the sur- 
face contribution 9 « Zq<j h = 2a(v t + Vb)/{\J~£t + ^/st) (x 
sgn(B)l&lQ. For S / the bulk contribution makes 
the dependence 9(B) nonmonotonic with the following 
asymptotics: 



Z a f 



1 + ZndZ 



< 1, 



Z a H , fir » 1, (19) 



The low-field limit (Or <C 1) agrees with the result of 
Ref. [13 which found 9 smaller than the surface con- 
tribution Zq<j h for nonzero bulk conductivity E, e.g. 
9 ss (T H /dY> <C Zq(T h for ZgdT, ^> 1. However, for strong 
fields VLt > 1 we find saturation of 9(B) precisely at the 
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FIG. 3: (a) Schematic geometry of EM transmission spec- 
troscopy of a TI film (see also text), (b) Magnetic-field 
dependent part of the transmission coefficient for a H 7^ 
(solid curves) and for a H = (dashed curves); lot — 1 and 
ZodT, = 1. (c) Same for lot = 35 and ZodE = 40. (d) Faraday 
angle in units of a = e 2 /ch versus magnetic field for zero and 
finite bulk conductivity E; lot <C 1. (e) Faraday angle versus 
frequency for zero and finite bulk conductivity E; Qt = 10. 
B and lo are in units of mc/er and 1/r, respectively. All data 
for e tj (, = 1 and thin film with ZodeLo/An ~ deui/c <g 1. 



surface value Zqu h because of the suppression of the bulk 
conductivity via classical cyclotron motion. This way of 
extracting the topological surface contribution may have 
an advantage over the previously proposed low-field de- 
tection scheme^ because the magnetization of the FMs in 
strong fields leads to a more robust surface Dirac gap A. 
As seen from Fig. the frequency dependence of the 
Faraday angle 9(iS) also saturates at the surface value 
Zq<t h , which can be used for its detection as well. 



VI. SUMMARY 

In summary, we have investigated galvanomagnetic 
and electromagnetic properties of topological insulators 
in which time-reversal symmetry is broken due to the sur- 
face quantum Hall effect. Our model includes both the 
dissipationless quantum Hall conductivity on the surface 
and the classical magnetoconductivity in the bulk of the 
system. Although the zero-field dc bulk conductivity may 
significantly exceed the surface one, the surface contri- 
bution can still be detected through anomalous magnetic 
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field dependencies of electrodynamic responses, reveal- 
ing the underlying broken time-reversal symmetry With 
appropriate modifications our findings can be extended 
to HgTe quantum wells which also support single- valley 
Dirac fermions^r— and show a pronounced Faraday ef- 
fect^ 
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